Plant Biotechnology JournalVolume 20, Issue 3 p. 420-422 Brief CommunicationOpen Access Profiling plant histone modification at single-cell resolution using snCUT&Tag Weizhi Ouyang, Weizhi Ouyang orcid.org/0000-0003-0099-4044 National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China These authors contributed equally to this article.Search for more papers by this authorShiping Luan, Shiping Luan National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China These authors contributed equally to this article.Search for more papers by this authorXu Xiang, Xu Xiang National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorMinrong Guo, Minrong Guo National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorYan Zhang, Yan Zhang National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorGuoliang Li, Guoliang Li orcid.org/0000-0003-1601-6640 National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China Hubei Key Laboratory of Agricultural Bioinformatics and Hubei Engineering Technology Research Center of Agricultural Big Data, 3D Genomics Research Center, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorXingwang Li, Corresponding Author Xingwang Li xingwangli@mail.hzau.edu.cn orcid.org/0000-0003-3985-2545 National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China *Correspondence (Tel +86-15007137931; fax +86-27-87384670; email xingwangli@mail.hzau.edu.cn)Search for more papers by this author Weizhi Ouyang, Weizhi Ouyang orcid.org/0000-0003-0099-4044 National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China These authors contributed equally to this article.Search for more papers by this authorShiping Luan, Shiping Luan National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China These authors contributed equally to this article.Search for more papers by this authorXu Xiang, Xu Xiang National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorMinrong Guo, Minrong Guo National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorYan Zhang, Yan Zhang National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorGuoliang Li, Guoliang Li orcid.org/0000-0003-1601-6640 National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China Hubei Key Laboratory of Agricultural Bioinformatics and Hubei Engineering Technology Research Center of Agricultural Big Data, 3D Genomics Research Center, Huazhong Agricultural University, Wuhan, ChinaSearch for more papers by this authorXingwang Li, Corresponding Author Xingwang Li xingwangli@mail.hzau.edu.cn orcid.org/0000-0003-3985-2545 National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China *Correspondence (Tel +86-15007137931; fax +86-27-87384670; email xingwangli@mail.hzau.edu.cn)Search for more papers by this author First published: 21 December 2021 https://doi.org/10.1111/pbi.13768Citations: 1 AboutSectionsPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onFacebookTwitterLinked InRedditWechat Chromatin immunoprecipitation with sequencing (ChIP-seq) with population cells or tissues yields ensemble epigenomic profiles that only represent the population average, which eliminates cell-to-cell epigenetic heterogeneity. Chromatin immunocleavage with sequencing (ChIC-seq) (Ku et al., 2019) that based on antibody-guided chromatin cleavage under targets is a practical alternative to ChIP-seq. Many single-cell ChIC-based methods, such as scChIC-seq (Ku et al., 2019), CUT&Tag (Kaya-Okur et al., 2019), CoBATCH (Wang et al., 2019), ACT-seq (Carter et al., 2019), and 10× scCUT&Tag (Bartosovic et al., 2021; Wu et al., 2021), have been developed and applied to study cell-type profiling, developmental trajectory, epigenetic heterogeneity, and transcriptional regulation in animals. However, due to the existence of cell walls, it is difficult to obtain single cells in plants. On the other hand, the existing scChIC methods require many single-cell barcoding procedures after tagmentation, which might lead to some DNA leakage and reduced mapping efficiency. Hence, a robust single-cell ChIC-seq method for plant epigenomic research is required emergently. In our previous study, we developed a rapid ChIC-based chromatin profiling protocol, nucleus CUT&Tag (nCUT&Tag) (Ouyang et al., 2021). Here, we combined the nCUT&Tag assay with 10× Single Cell ATAC, developing an easy-to-use single-nucleus CUT&Tag (snCUT&Tag) method in rice. We reformulated the bulk-cell nCUT&Tag with the droplet-based single-cell barcoding technology to develop snCUT&Tag (Figure 1a). Briefly, ~106 of fixed rice-seedling nuclei were collected and immuno-cleaved with protein G-Tn5 (Vazyme, Nanjing, China, cat no. TD901). Following Tn5-assisted in situ tagmentation, ~20 000 Tn5-indexed nuclei were loaded onto the 10× Genomics Single Cell ATAC microfluidics device, in which each nucleus is partitioned in a gel-bead-in-emulsion droplet and barcoded with unique oligonucleotides (Figure 1a). Finally, the chromatin features of every single nucleus could be captured by pooling and sequencing the barcoded DNA fragments. Figure 1Open in figure viewerPowerPoint snCUT&Tag profiles histone modification at single-cell resolution in rice. (a) Schematics of single-nucleus CUT&Tag workflow. (b) Representative H3K4me3 landscapes generated by indicated methods. Profiles of 3679 aggregate single nuclei (agg sc) and 80 selected single nuclei (80 sc) are presented. The genomic region contains 16 genes with gene IDs from OsMH_02G0056600 to OsMH_02G0058100. (c) Summary of the snCUT&Tag data. (d) Scatterplots showing Spearman correlations among snCUT&Tag, nCUT&Tag (nC&T), and eChIP-seq. (e) Violin plot showing FRiP scores of the snCUT&Tag data. Simulated random profiles are presented as control. (f) Heatmaps showing H3K4me3 peak signals across gene bodies. (g) Genomic screenshots showing aggregate marker peaks for 17 clusters (C1−C17). The genes IDs that annotated by marker peaks are presented in every row of the cell clusters. (h) Genomic screenshot showing the predicted enhancer–promoter interactions in indicated cell clusters and aggregate single cells (agg. predicted interactions). The green boxes indicate predicted enhancers (e1−e3). (i) Venn diagram showing overlapping of the predicted enhancers and the ATAC-seq peaks. (j–k) Density profiles and heatmaps showing enrichment of ATAC-seq signals around the predicted enhancers. We first profiled the H3K4me3 histone mark in seedlings from the Xian group rice (Minghui 63) using snCUT&Tag. We generated over 480 million read pairs, of which 98.7% were indexed with a unique barcode (Figure 1c). After mapping the reads to Minghui 63 reference genome and removal of duplicates, we obtained the H3K4me3 profiles of 3679 single nuclei, with a median of 401 fragments per nucleus (Figure 1c). For the 3679 aggregate single nuclei, we called 28 795 peaks (Figure 1b,c). Next, we compared the aggregate snCUT&Tag data with bulk-cell nCUT&Tag (Ouyang et al., 2021) and eChIP-seq data (Zhao et al., 2020). The aggregate snCUT&Tag peaks exhibited a high degree of consistency with nCUT&Tag and eChIP-seq data (Figure 1b). Coinciding with the high consistency, global scatterplots displayed significant correlations among the snCUT&Tag, nCUT&Tag, and eChIP-seq libraries (Figure 1d). As demonstrated by few reads falling in non-peak regions (Figure 1b), the snCUT&Tag data exhibited a significantly high fraction of reads in peaks (FRiP) value (median FRiP = 0.923) (Figure 1c,e). snCUT&Tag displayed similar peak signal profiles to nCUT&Tag and eChIP-seq, with H3K4me3 peaks mainly enriching at transcription start sites (Figure 1f). These results indicated that snCUT&Tag is a robust method that can be used to profile single-nucleus chromatin features. Since the epigenomic data can be used to predict chromatin interactions (Fulco et al., 2019; Ouyang et al., 2020), we examined the H3K4me3 snCUT&Tag data and predicted cell-type specific enhancer–promoter interactions according to the activity-by-contact model (Bartosovic et al., 2021; Fulco et al., 2019). Firstly, we performed dimensionality reduction and clustering according to Bartosovic et al. (2021). The results showed that rice seedlings were partitioned into 17 cell clusters, with significantly high enrichment-signal specificity (Figure 1g). Furthermore, we annotated cell types for the 17 clusters by investigating the marker peak-associated functional genes. For instance, a marker peak in cluster 6 showed specific H3K4me3 enrichment around a malate dehydrogenase-coding gene, which was reported as a marker gene of mesophyll cells in maize (Marand et al., 2021), implying that cluster 6 should be classified as mesophyll cells. Then, we predicted enhancer–promoter interactions. Our results demonstrated that the predicted chromatin interactions were well-overlapped with H3K4me3 ChIA-PET data, and many enhancer-centred chromatin interactions showed cell-type biases that only occurred in specific cell clusters (Figure 1h). Expectedly, ~51.5% (15 881 of 30 839) of predicted enhancers overlapped with ATAC-seq peaks (Figure 1i), showing significant enrichment of ATAC-seq signals around the centre of the predicted enhancers (Figure 1j,k). Our results indicate that snCUT&Tag can be applied to dissect epigenetic heterogeneity and predict cell type-specific enhancer–promoter interactions, which might facilitate the annotation of regulatory elements and reconstruction of 3D genome structures. In conclusion, we developed a simple and practical single-nucleus ChIC-seq method (snCUT&Tag) for single-cell plant epigenomic studies. We demonstrated that snCUT&Tag could be used to dissect the cell lineage and epigenetic heterogeneity and reconstruct chromatin topology of each single nucleus. Since nCUT&Tag is applicable for studying bulk-cell chromatin features with nuclei isolated from various plant species and tissues (Ouyang et al., 2021), we believe that snCUT&Tag may have a broad spectrum of application in profiling single-nucleus histone marks for other types of plant materials as well. Expectedly, snCUT&Tag provides an avenue for plant single-cell epigenomic studies and may help understand the 3D genomic and epigenetic basis of transcription regulation in a single cell. Acknowledgements This work was supported by the National Natural Science Foundation of China (32070612, 31771422), the funds from the Interdisciplinary Sciences Research Institute (2662021JC005), and the National Key Laboratory of Crop Genetic Improvement (ZK201906), as well as the Baichuan Project (to W.O.) from College of Life Science and Technology, Huazhong Agricultural University. Conflict of interest No conflict of interest declared. Author contributions X.L. and G.L. supervised the research. W.O. performed the experiments with assistance from X.X., M.G., and Y.Z. S.L. and W.O. performed bioinformatic analysis. W.O. and X.L. wrote the manuscript. Open Research Data Availability Statement The snCUT&Tag data have been deposited in the National Genomics Data Center (https://ngdc.cncb.ac.cn/) with accession number CRA004386. References Bartosovic, M., Kabbe, M. and Castelo-Branco, G. (2021) Single-cell CUT&Tag profiles histone modifications and transcription factors in complex tissues. Nat. Biotechnol. 39, 825– 835. CrossrefCASPubMedWeb of Science®Google Scholar Carter, B., Ku, W.L., Kang, J.Y., Hu, G., Perrie, J., Tang, Q. and Zhao, K. (2019) Mapping histone modifications in low cell number and single cells using antibody-guided chromatin tagmentation (ACT-seq). Nat. Commun. 10, 3747. CrossrefPubMedWeb of Science®Google Scholar Fulco, C.P., Nasser, J., Jones, T.R., Munson, G., Bergman, D.T., Subramanian, V., Grossman, S.R. et al. (2019) Activity-by-contact model of enhancer-promoter regulation from thousands of CRISPR perturbations. Nat. Genet. 51, 1664– 1669. CrossrefCASPubMedWeb of Science®Google Scholar Kaya-Okur, H.S., Wu, S.J., Codomo, C.A., Pledger, E.S., Bryson, T.D., Henikoff, J.G., Ahmad, K. et al. (2019) CUT&Tag for efficient epigenomic profiling of small samples and single cells. Nat. Commun. 10, 1930. CrossrefPubMedWeb of Science®Google Scholar Ku, W.L., Nakamura, K., Gao, W., Cui, K., Hu, G., Tang, Q., Ni, B. et al. (2019) Single-cell chromatin immunocleavage sequencing (scChIC-seq) to profile histone modification. Nat. Methods, 16, 323– 325. CrossrefCASPubMedWeb of Science®Google Scholar Marand, A.P., Chen, Z.L., Gallavotti, A. and Schmitz, R.J. (2021) A cis-regulatory atlas in maize at single-cell resolution. Cell, 184, 3041– 3055. CrossrefCASPubMedWeb of Science®Google Scholar Ouyang, W., Xiong, D., Li, G. and Li, X. (2020) Unraveling the 3D genome architecture in plants: present and future. Mol. Plant, 13, 1676– 1693. CrossrefCASPubMedWeb of Science®Google Scholar Ouyang, W., Zhang, X., Peng, Y., Zhang, Q., Cao, Z., Li, G. and Li, X. (2021) Rapid and low-input profiling of histone marks in plants using nucleus CUT&Tag. Front. Plant Sci. 12, 634679. CrossrefPubMedWeb of Science®Google Scholar Wang, Q., Xiong, H., Ai, S., Yu, X., Liu, Y., Zhang, J. and He, A. (2019) CoBATCH for high-throughput single-cell epigenomic profiling. Mol. Cell, 76, 206– 216.e7. CrossrefCASPubMedWeb of Science®Google Scholar Wu, S.J., Furlan, S.N., Mihalas, A.B., Kaya-Okur, H.S., Feroze, A.H., Emerson, S.N., Zheng, Y. et al. (2021) Single-cell CUT&Tag analysis of chromatin modifications in differentiation and tumor progression. Nat. Biotechnol. 39, 819– 824. CrossrefCASPubMedWeb of Science®Google Scholar Zhao, L., Xie, L., Zhang, Q., Ouyang, W., Deng, L., Guan, P., Ma, M. et al. (2020) Integrative analysis of reference epigenomes in 20 rice varieties. Nat. Commun. 11, 2658. CrossrefCASPubMedWeb of Science®Google Scholar Volume20, Issue3March 2022Pages 420-422 FiguresReferencesRelatedInformation